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Abstract 

We introduce a new smearing transformation, the hypercubic (HYP) fat link. The 
hypercubic fat link mixes gauge links within hypercubes attached to the original link only. 

Using quenched lattices at /? = 5.7 and 6.0 we show that HYP fat links improve flavor 
symmetry by an order of magnitude relative to the thin link staggered action. The static 

potential measured on HYP smeared lattices agrees with the thin link potential at 
distances r/a > 2 and has greatly reduced statistical errors. These quenched results will be 
used in forthcoming dynamical simulations of HYP staggered fermions. 
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I. INTRODUCTION 



Fat link actions, actions where the thin link fermionic gauge connection is replaced by 
a combination of extended but local gauge paths in a gauge invariant manner, have been 
used in a variety of quenched and dynamical simulations. Since the gauge paths of fat 
link actions combine only a finite number of gauge links in the local environment of the 
original link, fat link actions are in the same universality class as their thin link counterpart. 
For many physical observables fat link actions show reduced scaling violations. Staggered 
fat link fermions have improved flavor symmetry, both in quenched ||1], 0, ^ and dynamical 
simulations Wilson-type clover fat link fermions have been studied only in quenched 

simulations where they show improved chiral behavior resulting in better scaling of the 
hadron spectrum and fewer exceptional configurations |0, If. Fat links have been used in 
overlap simulations as a starting action for evaluating Neuberger's formula |||, |l^, 0. The 
improved chiral properties of fat links result over an order of magnitude improvement in 
convergence. 

It is quite intuitive to see why fat link staggered fermions show improved flavor behav- 
ior. Flavor symmetry breaking of staggered fermions is the consequence of distributing the 
components of the Dirac spinors to different lattice sites within a hypercube. Since each 
lattice site couples to different gauge fields, the different flavor and Dirac components feel 
different gauge environments. The original SU{nf) x SU{nf) symmetry is explicitly broken 
to a remnant U{1) symmetry. Unifying the gauge fields within the hypercube would solve 
this problem. Smearing, while does not unify the gauge fields, can remove some of the local 
fluctuations thus improving flavor symmetry. The most effective smearing is one that max- 
imally smoothes the gauge fields within the hypercube but does not extend much beyond 
it. 

The situation is similar for Wilson fermions that break chiral symmetry explicitly. The 
Wilson Dirac operator is not anti-Hermitian, its eigenvalues are not constrained to the imag- 
inary axis. The real eigenvalues spread around Kcr causing spurious poles in the quenched 
quark propagator in the positive pion mass region |0, ||, Since the real eigenvalue 

modes of the Dirac operator are instanton related, one can study the spread of the real 
eigenvalues on smooth instanton backgrounds. The real eigenmodes of the Wilson Dirac 
operator move farther from K^r into the positive mass region on small instantons 0. The 
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Dirac operator has real modes even on configurations where the "instanton", obtained by re- 
peated renormalization group blocking from a large classical solution, is small with respect 
to the lattice spacing. These objects do not have topological interpretation, they are only 
dislocations, yet they are responsible for the worst spurious poles of the Dirac propagator. 
On ordinary Monte Carlo gauge configurations dislocations arise from large fiuctuations of 
the plaquette. Just as dislocations distort the topological charge measurements they also 
distort the low lying Dirac spectrum and any physical observable related to it. Smearing 
removes some of the fluctuations and their corresponding spurious eigenmodes. This reduces 
the occurence of exceptional configurations, improves scaling [0, ||] and makes topological 



charge measurement on pure gauge configurations possible ||TJ, [T^. Again, the most effective 
smearing is one that maximally smoothes the gauge fields within a small, hypercube sized 
volume, but does not extend much beyond it. 

There are observables where some fat link actions have unacceptably large scaling vio- 
lations. Fattening changes the local, short distance properties of the gauge configurations. 
On lattices with large lattice spacing too much fattening can destroy short distance physical 
properties. This is most evident in the Coulomb sector of the static potential ll^. Increased 
levels of smearing destroy more and more of the short distance potential though it does 
not change the long distance string tension. The heavy-light decay constants are also very 
sensitive to smearing. Even low levels of smearing of the heavy fermions change the value 
of /b significantly ||T^, |T^. Heavy fermions feel the short distance structure of the vacuum 
and that gets distorted by the fattening. If only the light quark action is based on fat links, 
this problem is avoided while the chiral properties of the light quarks are improved. These 
expectations are supported by preliminary results The fact that certain observables 

show large scaling violations with fat link actions does not mean that fat link actions are 
incorrect, they are just not effective in every situation. 

It would be desirable to construct fat links that are local enough not to destroy short 
distance quantities yet smooth enough to show the desired improved chiral or flavor behavior 
of fat link actions. In this paper we describe a smearing transformation, hypercubic blocking, 
along these lines. This transformation improves flavor symmetry as much as three levels of 
APE blocking but mixes links from hypercubes attached to the original link only. The static 
potential changes only at small distances r/a < 2 where lattice artifacts are anyhow large. 
This paper deals with the definition and quenched properties of the hypercubic fat link 
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action. In a recent paper we showed how to simulate fat link dynamical systems. We 
were able to simplify and modify this algorithm in such a way that dynamical simulations 
with the hypercubic action are possible and efficient. We used the dynamical hypercubic 
action to study finite temperature QCD with four flavors of staggered quarks |I9| . 

This paper is organized as follows. In Sect. 2 we define the hypercubic blocking. The 
transformation has three tunable parameters. We describe how the parameters are optimized 
and show that as the lattice spacing decreases the transformation gets more effective, i.e. the 
number of smearing levels can be kept fixed and small even as one approaches the continuum 
limit. We show that the hypercubic action considerably improves the flavor symmetry of 
staggered quarks in quenched simulations. The relative improvement with respect to the 
standard action is about an order of magnitude and increases as the lattice spacing decreases. 
In Sect. 3 we present quenched results for the static potential. We show that the hypercubic 
action allows the determination of the Sommer scale rg and the string tension a with better 
accuracy compared to calculations done with the thin link standard action. The distortion 
of the potential at small distances due to the hypercubic blocking is limited to distances 
r/a < 2 independent of the lattice spacing and rotational symmetry is also improved. In 
Sect. 4 we summarize our results. 



II. THE HYPERCUBIC BLOCKING 



A. Definition and Optimization of the Hypercubic blocking 

The fat links of the hypercubic blocking (HYP) are constructed in three steps. At the 
final level the blocked link I^,/^ is constructed via projected APE blocking [^] from a set of 
decorated links Vi^^j_-y as 

where Ui^^ is the original thin link and the index z/ in Vi,^;:^ indicates that the fat link 
at location i and direction is not decorated with staples extending in direction z/. The 
decorated links Vi^^-^ are constructed with a modified projected APE blocking from an other 
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a) 



b) 




FIG. 1: Schematic representation of the hypercubic blocking in three dimensions, a) The fat hnk is 
built from the four double-lined staples, b) Each of the double-lined links is built from two staples 
which extend only in the hypercubes attached to the original link. An important point is that the 
links entering the staples are projected onto SU(3). 



set of decorated links, as 



(2) 



where the indices pu indicate that the fat link Vi^^-pi, in direction /i is not decorated with 
staples extending in the p or z/ directions. The decorated links Vi^^.^py are constructed from 
the original thin links with a modified projected APE blocking step 



±v¥=p,i^<p 



rt 



(3) 



Here only the two staples orthogonal to /i, u and p are used. With the construction eqs. 

the fat link Vi^p mixes thin links only from hypercubes attached to the original link. 
The hypercubic blocking is schematically represented in figure |1|. The parameters ai, 02 and 
^3 can be optimized to achieve the smoothest blocked link configuration. The construction 
eqs. (|l])-(||) can be iterated. 

We could use a physical quantity like the level of fiavor symmetry restoration to optimize 
the parameters of the blocking, but in a three dimensional parameter space that approach 
is very computer intensive. Instead we have decided to look at the underlying reason of 
fiavor and chiral symmetry violation. As we explained in the introduction, both the stag- 
gered fermions' fiavor symmetry violation and Wilson fermions' exceptional configurations 
can be understood as fiuctuations in the local plaquette. It is not the average, not even the 



50 




a) 



JTruwrillflTflj 1 1 1 1 1 1 ii 



-1 1 

Min(DU) 




-1 1 
Mm(DU) 



FIG. 2: a) The histogram of the smallest plaquette distribution on a set of 4^ /? = 5.7 configurations 
(dashed lines) and after hypercubic blocking (solid lines) on the same configurations, b) Same as 
a) but the solid lines correspond to one level APE smearing with a = 0.7. 

average fluctuation of the plaquette that causes the problems but the few largest plaquette 
fluctuations that create dislocations. To minimize the largest fluctuations we have to max- 
imize the smallest plaquettes of the system. We do this by considering the distribution of 
the smallest plaquette in a finite, fixed volume lattice. The optimized hypercubic blocking 
parameters will give the largest expectation value for the smallest plaquette in this finite 
volume ensemble. We have optimized the a parameters using a set of 500 4^* lattices created 
at P = 5.7 with the Wilson plaquette action. 

Figure |2|/a shows the histogram of the smallest plaquette on each of the 500 original 
configurations (dashed lines) and after hypercubic blocking (solid lines) with parameters 



ai = 0.75, 0^2 = 0.6, = 0.3. 



(4) 



The average of the smallest plaquette jumps from -0.75 to 1.95. This is about the best 
that can be achieved with the hypercubic blocking. It is important to remember that we 
are maximizing the smallest plaquette values on the lattice, not the average plaquettes, 
which on these lattices are 1.67 and 2.82, respectively. The minimum plaquette is not 
a well defined physical quantity, it strongly depends on the lattice size, but as long as 
we compare different actions on the same set of configurations, the minimum plaquette 
shows how the most severe vacuum fiuctuations change during smoothing. In figure |2|/b 
we compare the minimum plaquette distribution of the thin link action to a one-level APE 
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FIG. 3: The parameter 82 of eq. (P) for the two hghtest non-Goldstone pions VTj 5 (lower points) and 
TTjj- (higher points) as the function of the minimum plaquette average of the fat hnk actions. The 
pion spectrum is computed on 8^ x 24, (3 = 5.7 quenched configurations, the minimum plaquette 
average on 4^, /? = 5.7 configurations. The dotted (vTj^s) and dashed (vTjj) lines show how the flavor 
symmetric limit is approached as the minimum plaquette approaches three. 

smeared action with a = 0.7 on the same configurations that were used to obtain figure 
H/a. Here the minimum plaquette average increases only to 0.45 after smearing. Note that 
it is very important to project to SU(3) the intermediate decorated links in eqs. (|l])-(0)- 
Without the projection the resulting hypercubic smearing is not better than a one-level APE 
transformation. 



B. Flavor symmetry Restoration with Fat Link Actions 

In the previous section we described an intuitive and fast way to optimize fat link actions 
by maximizing the minimum plaquette average in finite volume lattices. In this section we 
show that, as expected, this method optimizes fiavor symmetry as well. We compare the 
level of fiavor symmetry restoration of a few fat link actions with their respective minimum 
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plaquette average computed on 4"^ configurations at /3 = 5.7. We use the parameter ||^ 

(5) 



ni?p — tuq 



extrapolated to mG/rrip = 0.55 to quantify the diff'erence between the mass of the true 
Goldstone pion mo and the mass of the pseudoscalar particles or pions m^r. For a flavor 
symmetric action ^2 = for all the pions vr. We calculated the meson spectrum on an 
ensemble of 8^ x 24, = 5.7 quenched configurations using the thin link, the N=l, 2, and 
3 levels probabilistic APE blocked actions with a = 0.7 and projection parameter A = 500 
(for the definition of probabilistic blocking and the projection parameter see Ref. |^), the 
N=l deterministic (A = oo) APE blocked action and the hypercubic blocked action. Our 
results are summarized in figure | where we plot 62 for the two lightest non-Goldstone 
pions, TTj 5 with flavor structure 7^75 (lower points) and iTij with flavor structure 7j7j (higher 
points) as the function of the minimum plaquette of the fat link actions. 62 seems to depend 
on the minimum plaquette linearly approaching the flavor symmetric limit ^2 = as the 
dotted (vTj 5) and dashed (vTj j) lines in figure | indicate, though we do not have a theoretical 
explanation why this dependence should be linear. The HYP action has about the same 
flavor symmetry violation than the N=3 APE action that we studied in Ref. |||] reducing S2 
by about a factor of six relative to the thin link action. 

The data in flgure |^ are extrapolated to mc/mp = 0.55, a value considerably higher than 
the physical value 0.175 of this ratio. Since flavor symmetry violation increases as the the 
quark mass decreases, it is important to see how far one can lower the quark mass before 
flavor symmetry is ruined. In flgure |[/a we plot the relative mass splitting 

A. = (6) 

rriG 

between the lightest non-Goldstone pion tTj 5 and the Goldstone pion G as the function of 
(romc)^ on /? = 5.7 quenched conflgurations both for the thin and HYP actions. The flrst 



three points at the lowest romc values for the thin link action are from Ref. IgT]. The 
ratio diverges as the quark mass approaches zero but this divergence is much slower for 
the HYP action. The relative mass splitting for the HYP action grows to about 30% at 
fofnc = 0.70 where the ratio mc/mp of the Goldstone pion to the rho mass is 0.23(8). At 
similar values of the mo j rap ratio the mass splitting with the thin link action is between 
130% and 190%. 
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FIG. 4: a) The relative mass splitting A,r of eq. @ between the lightest non-Goldstone pion tTj^s 
and the Goldstone pion G as the function of (romc)^ on /? = 5.7 quenched configurations both for 
the thin (fancy diamonds) and the HYP (bursts) actions, b) Same as a) but on (3 = 6.0 quenched 
configurations where the lattice spacing is reduced by about a factor of 2. The first three points at 



the lowest romc values for the thin link action in both figures are from Ref. |21|. 



C. Hypercubic blocking on Finer Lattices 

We have demonstrated that the minimum plaquette value strongly correlates with the 
level of flavor symmetry restoration at the same lattice spacing. The question naturally 
arises: if we decrease the lattice spacing, will the smearing be effective or do we have to 
increase the number of smearing iterations? In figure |^ we show the minimum plaquette 
value calculated with the thin and hypercubic actions with the optimized parameters of eqn. 
(I) at P = 5.7,6.0 and 6.2. These calculations were done on 8^ lattices, that is why the 
P = 5.7 result looks significantly worse than before: on a larger lattice the probability of 
finding a large plaquette fluctuation is higher. The difference between the three /? values is 
striking. One certainly does not need more smearing at a higher coupling. As /3 increases, the 
fat link minimum plaquette value increases faster than the corresponding thin link minimum 
plaquette. That implies that one level of hypercubic blocking will be effective for all coupling 
values. To verify this expectation we calculated the HYP staggered meson spectrum on a set 
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FIG. 5: The histograms of the smallest plaquette distribution on sets of 8^ gauge configurations at 
different (3 values (dashed lines) and after hypercubic blocking (solid lines) on the same configura- 
tions. 

of 16^ X 32, /3 = 6.0 configurations. The configurations were downloaded from the NERSC 



archive |^. The mass splitting for the lightest non-Goldstone pion vTj 5 is plotted as the 
function of (rom^)^ in figure ||/b. The thin link data are from Ref. The vertical scale 
of figure §/b is one forth of the vertical scale of ||/ a but the qualitative features of the two 
plots are very similar. decreases by an order of magnitude from thin to HYP action, 
fairly independent of the lattice spacing or quark masses. 



III. THE STATIC POTENTIAL OF FAT LINK ACTIONS 



In the following we consider the static potential, an observable that is quite sensitive to 
the smearing of the gauge links. We analyzed two sets of lattices. The first set consists 240 
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8^ X 24 configurations generated with the Wilson plaquette action at /3 = 5.7. The second 
set is an ensemble of 220 16'^ x 32 lattices from the NERSC archive generated also with 



Wilson plaquette action at /3 = 6.0 [^]. We calculate the static quark potential on both 
sets using thin links and hypercubic smeared links. 

The static potential on pure gauge lattices has been measured and analyzed by several 
groups on larger lattices with orders of magnitude higher statistics. Our goal here is not 
the improvement of existing potential measurements, rather we want to demonstrate that 
the HYP action changes the potential only at small distances, r/a < 2 and, due to the 
smoothness of the gauge configuration, provides an improved signal both for the Sommer 



parameter tq |^ and even more for the string tension. The HYP potential also shows 
reduced rotational symmetry violation. Our fitted values for r^/a and ay/a are consistent 
with previously published results. 

We extracted the potential by fitting the time dependence of the different Wilson loops to 
a single exponential. We have measured the on-axis (m,0,0), m < n.j;/2 and three different 
off-axis sets: (m,m,0), m < n^j^^ (m,m,m), m < n^j^ and (m,2m,0), m < 12^/4: {rix is 
the number of lattice points in the spatial directions). To reduce statistical errors we used 
several (20-30) levels of transversal APE smearing on the thin link lattices. After hyper- 
cubic blocking the configurations are so smooth that transversal smearing is not important 
anymore. We used only 5-10 levels of transversal smearing. The statistical errors are from 
standard jackknife analysis. 

In figures ^ and |^ we plot both the thin link and HYP potentials for the /5 = 5.7 and 
[3 = 6.0 data sets. The thin link potentials are shifted by a constant value to match the HYP 
potentials in the 2.5 < r/a < 5.5 (/? = 5.7 ) and 2.5 < r/a < 8 (/? = 6.0) range. Even though 
the first 5 r/a values were not used for the matching in either coupling, the potentials match 
closely for r/a > -\/2. The distortion effects in the potential due to smearing are observable 
only in the first two r/a values, independent of the lattice spacing. Moreover, the HYP 
potentials have significantly reduced statistical errors, about a factor of 4-7 at /? = 5.7, 
10-15 at P = 6.0. This is most obvious at large distances. For example the /? = 6.0 thin 
link potential signal becomes extremely noisy for r/a > 7 with our statistics while the HYP 
potential can be easily measured even at the largest distance, r/a = 11.31 . 
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FIG. 6: The thin hnk and HYP potentials for the /? = 5.7 data set. The thin hnk potential is shifted 

by a constant value to match the HYP potential in the 2.5 < r/a < 5.5 range. The solid line is 
identical in both plots. It is a three parameter fit of the HYP potential in the range r/a G (1.9, 5.5). 
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FIG. 7: The thin link and HYP potentials for the (3 = 6.0 data set. The thin link potential is 
shifted by a constant value to match the HYP potential in the 2.5 < r/a < 8 range. The solid line 
is identical in both plots. It is a three parameter fit of the HYP potential in the range r/a G (2.5, 12). 
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TABLE I: Results obtained with the HYP potential using a global fit of the form of eq. (0). Also 
included are some large scale simulation results at the same coupling values. 



To extract ro and y/a from our data we fit the potential in the form 

Co 



V{r) 



+ Ci + C2r 



(7) 



in a range r^m < r < r^ax- fmax is limited by the accuracy of our data, Vmax/O' = 6 — 7 
for the thin link, 10 — 12 for the HYP potential at /? = 6.0. We choose rmin to exclude 
the region where the potential is changed by the smearing or where short distance lattice 
artifacts are large, typically rmm/a = 1 — 3. We do not include any lattice correction in 
the fit. For one, we do not know the perturbative potential for the HYP action, though 
it can be calculated. On the other hand we do not seem to need it as the HYP potential 
data shows only minimal lattice artifacts. Our results obtained with the HYP potentials 
are summarized in Table |I[ Fits to our thin link potentials gave consistent results but with 
much higher statistical errors and are not listed. In Table | we also include results from a 



few large scale potential simulations with thin links ||2J, |2^. We considered several Vmax and 
rmin values and included all on- and off-axis data points in the given range in a global fit to 
eq. (|^). The errors in Table | represent the average jackknife error of the different fits. Our 
values for rg/a and ay/a from the HYP potential are consistent with the results from large 
scale computations. 

We would like to emphasize that the HYP potential has about an order of magnitude 
smaller statistical errors than the thin link potential on the same set of configurations. This 
allowed us to extract the parameters of the potential with precision comparable to the large 
scale simulations even though our analysis is not as refined. In particular we did not use 



the one-link integral or multi-hit procedure ||26|| for the timelike links in the Wilson loops as 
has been done in ||2^ . 
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FIG. 8: The HYP potential at /3 = 6.0. The dashed lines are fits with eq. (|7|) along the different 
on- and off-axis directions the potential is measured. All four fits use the range r/a G (2.5, 12). 

In addition to the reduced statistical errors the HYP potential shows better rotational 
invariance as well. We use the same quenched configurations both for the thin and HYP 
potentials, but rotational symmetry violations can come not only from the configurations 
but also from the operator used to measure the potential. As an example we mention the 
potential for a perfect action [2^. A perfect action has no scaling violations, yet the pertur- 
batively calculated potential is not rotational invariant unless one uses a perfect potential 
operator. The HYP operator is of course not perfect, but it appears to be better than the 
thin link Wilson loop operator. To illustrate rotational symmetry restoration, we separated 
the on-axis and the three off-axis directions of the (3 = 6.0 potential and fitted them inde- 
pendently with the form of eq. (|^). The fits for the thin link potential were very poor but 
for the HYP potential all four fits gave consistent results. We fitted the HYP potential in 
the range r/a E (2.5, 12). The fitted values for ro/a were all within one standard deviation, 
for ay/a within two standard deviations. The fit results are shown in figure ^ The four 
fitted lines agree very well in the range r/a G (2.5,6.5), they deviate at larger r/a due to 
the slightly different ay/a values. The uncertainty of the fitted lines is comparable to their 
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spread in figure |^. Closer inspection of the data points reveals some rotational symmetry 
violation even with the HYP potential. It should be possible to correct for some of these 
effects perturbatively thereby increasing the accuracy of the scale determination. 

The improved rotational symmetry of the HYP potential is not the main motivation for 
using the HYP operator to measure the potential. Improved gauge actions can be even 
more effective in restoring rotational invariance of the static potential |2^, The 
improvement of statistical accuracy by an order of magnitude of the HYP potential as 
compared to the thin link potential is a more important feature. 

IV. SUMMARY 

In this paper we introduced a new kind of fat link, the hypercubic (HYP) fat link. The 
HYP fat link is constructed to be localized, it does not extend beyond hypercubes attached 
to the original link. The construction has three free parameters that are optimized to reduce 
large plaquette fluctuations. 

Using quenched conflgurations at /3 = 5.7 and (3 = 6.0 we calculated the staggered 
spectrum of the HYP action. We showed that at both coupling values flavor symmetry is 
greatly improved. The relative mass splitting, = {uig — mT,)/mG is reduced by an order 
of magnitude at flxed tquig-, independent of the coupling and the quark mass as well. 

We have also considered the static potential calculated on HYP smeared /5 = 5.7 and 
(3 = 6.0 conflgurations. We showed that smearing changes the potential only at distances 
r/a < 2 yet reduces statistical errors by about an order of magnitude. Even with the very 
limited statistics of 220-240 lattices we could extract both the Sommer parameter r^/a and 
the string tension ay/a. Our results for both of these quantities are consistent with large 
scale simulation results with errors that are better than our statistics would indicate from 
the thin link potential. 

We have undertaken this quenched study on the properties of the hypercubic blocking 
because we want to do dynamical simulations with the HYP action. Such simulations are 
feasible using combined over-relaxation and heat-bath updates. In a forthcoming publication 
ITqII we use the dynamical HYP action in a flnite temperature simulation to show how 
improved flavor symmetry changes the phase diagram of four-flavor dynamical staggered 
fermions. 
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